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ABSTRACT 

In giant molecular clouds (GMCs), shocks driven by converging turbulent flows 
create high-density, strongly-magnetized regions that are locally sheetlike. In previous 
work, we showed that within these layers, dense filaments and embedded self-gravitating 
cores form by gathering material along the magnetic field lines. Here, we extend the 
parameter space of our three-dimensional, turbulent MHD core formation simulations. 
We confirm the anisotropic core formation model we previously proposed, and quantify 
the dependence of median core properties on the pre-shock inflow velocity and upstream 
magnetic field strength. Our results suggest that bound core properties are set by the 
total dynamic pressure (dominated by large-scale turbulence) and thermal sound speed 
Cg in GMCs, independent of magnetic field strength. For models with Mach number 
between 5 and 20, the median core masses and radii are comparable to the critical 
Bonnor-Ebert mass and radius defined using the dynamic pressure for Pgxt- Our results 
correspond to Mcore = si^nd i?core = 0.34c5^(Gpo'^o^)~^^^ po si^nd 

vq the large-scale mean density and velocity. For our parameter range, the median 
Mcore ^ 0.1 — IMq, but a very high pressure cloud could have lower characteristic core 
mass. We find cores and filaments form simultaneously, and filament column densities 
are a factor ^ 2 greater than the surrounding cloud when cores first collapse. We also 
show that cores identified in our simulations have physical properties comparable to 
those observed in the Perseus cloud. Superthermal cores in our models are generally 
also magnetically supercritical, suggesting that the same may be true in observed clouds. 

Subject headings: ISM: magnetic fields — MHD — turbulence — stars: formation 


Introduction 


Prestellar core formation in giant molecular clouds (GMCs) is an important issue in theoretical 


studies of star formation, because these cores are the immediate precursors of protostars (Shu et 


al. 1987; McKee & Ostriker 2007; Andre et al. 2014). It is believed that the magnetic field and 
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supersonic turbulence may both play important roles in core formation and subsequent evolution. 
In GMCs, simulations starting more than a decade ago have shown that overdense structures 
generated by supersonic turbulence may collapse gravitationally to form protostellar systems, while 


also attracting material from their surroundings (e.g. 

Ballesteros-Paredes et al.||1999; Ostriker et al. 

1999 

; Klessen et al.|2000; 

Padoan et al.|2001 

; Bate et al.|2003). Magnetic fields limit compression in 


large-scale turbulence-induced shocks, channel material toward forming filaments, provide support 


for cores as they grow, and remove angular momentum in collapsing cores ( 

Mestel & Spitzer||1956; 

Strittmatter||1966HMouschovias & Spitzer |1976l Mestel||1985; |Mouschovias||1991 ; 

Allen et al.||2003l 

Basu et al.||2009; 

Li et al.||2010, 

2014). 


Because GMCs are only lightly ionized, and magnetic fields are only coupled to charged par¬ 
ticles, magnetic stresses are mediated by ion-neutral collisions, and are affected by the level of am- 
bipolar diffusion. Analytic studies and numerical simulations have shown that supersonic motions 
accelerate ambipolar diffusion ( Fatuzzo fc Adam"s||2002 ; Heitsch et al.||2004 ; Li fc Nakamura|[2004 ). 
Similar simulations with both strong turbulence and ambipolar diffusion have also demonstrated 


core evolution times, efficiency of star formation, and core structure similar to observations (Naka¬ 


mura fc Li|2005 , 2008; Kudoh fc Basu|20^ , 2011; Basu et al.|2009 ). More recently, Chen & Ostriker 


(2012) studied the one-dimensional C-type shocks and identified a transient stage of turbulence- 
accelerated ambipolar diffusion. This transient stage, with timescale ^transient "^0.1 — 1 Myr (de¬ 
pending on ionization), can explain the enhanced diffusion rate and affect the magnetization of 
cores that form. 


In Chen &; Ostriker (2014, hereafter COM), we applied three-dimensional numerical simula¬ 
tions to study the roles of magnetic fields and ambipolar diffusion during prestellar core formation in 
turbulent cloud environments. Our simulations adopted the model framework of |Gong fc Ostriker 


(2011) to focus on the shocked layer produced by turbulent converging flows, and surveyed vary¬ 


ing ionization and angle between the upstream flow and magnetic field. In simulations, we found 
hundreds of self-gravitating cores with masses M ^ 0.04 — 2.5 Mq and sizes L ^ 0.015 — 0.07 pc, 
all formed within 1 Myr. 


In |C014| we also found that core masses and sizes do not depend on either the ionization or 
upstream magnetic field direction, and ambipolar diffusion is in fact not necessary to form low-mass 
supercritical cores. Our analysis showed that this is the result of anisotropic contraction along field 
lines, which can be clearly seen in our simulations, with or without ambipolar diffusion (see also 
Van Loo et al.|[Ml4 ). In the anisotropic core formation model, low-mass magnetically supercritical 
cores form rapidly even in a strongly magnetized medium with high ionization. This explains the 
prevalence of magnetically supercritical cores in observations ( Crutcher||2012 ). 

Using a simple scaling argument, |C014| suggested the characteristic core mass may be set by 
the mean turbulent pressure in a CMC, regardless of magnetic effects. The predicted core mass 
is a factor ^ lower than the thermal Bonnor-Ebert mass at the mean density in the cloud. 


where A4 is the turbulent Mach number. Previously, Gong & Ostriker (2011) proposed a similar 
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formula based on the preferred scale for gravitational fragmentation of post-shock layers, for the 


purely hydrodynamic case. Padoan et al. (1997) also argued for a similar characteristic mass, based 


on statistics of turbulent flows. Although the analyses of Gong & Ostriker (2011) and Padoan et al. 


(1997) neglect magnetic fields, the end result is similar to the prediction of COM that incorporates 


magnetic effects, with the turbulent pressure in a cloud setting the characteristic core mass. 


Here, following |C014 , we continue our study of anisotropic core formation in turbulent molec¬ 
ular clouds. We extend our previous numerical study to explore how the turbulent and magnetic 
pressures of the pre-shock gas can affect core formation in the compressed region. We demonstrate 
that the dependence of core properties on pre-shock parameters are similar to those predicted by 
the anisotropic core formation model of |C014 We also compare our results with observations. 


showing that the mass-size relationship and ratio of mass to critical value of our simulations is 


comparable to that seen in Perseus and other star-forming regions (Sadavoy et al. 2010a; Kirk et 


al. 2013). 


The outline of this paper is as follows. We review the anisotropic core formation model in Sec¬ 
tion outlining the successive dynamical stages and associated parameter dependence expected. 
Section describes the equations solved in our numerical simulations, and specifies the model pa¬ 
rameter set we shall consider. The post-shock gas structure for varying parameters is analyzed in 
Section]^ including physical properties of the compressed layer (Section |4.1[ ), and development of 
filaments within it (Section [4.2[ ). In Section]^ we provide quantitative results for masses, sizes, mag¬ 
netizations, and other physical properties of the bound cores identified from our simulations, and 
compare to predictions from COM We also compare these results with observations (Section]^, 
focusing on interpreting the physical state of super-Jeans mass cores and mass-size relationships. 
Section [3 summarizes our conclusions. 


2. Anisotropic Core Formation: Review 


Here we briefly review the anisotropic condensation model of core formation proposed in COM 


We consider a strongly-magnetized post-shock region created by a large-scale converging turbulent 
flow within a cloud. As shown in COM| (see Figures 3—5 there), the magnetic field will lie primarily 
parallel to the shock front in the layer no matter what the initial inclination angle is, because only 
the component of the field parallel to the shock front (or perpendicular to the inflow) is amplified. 
Equation (4) of C014 gives the compression ratio for the component of the magnetic field that 
is amplified. The anisotropic condensation model describes core formation in the post-shock layer 
with density pps and threaded by magnetic field strength Bps (Figuretop left). For a cylinder 
with radius R and length L along the magnetic field (Figure top right) ^ if 2R ^ L < Tmag,crit 
(see Equations (30) and (31) of COM) for 


Pps 27V 


( 1 ) 
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Fig. 1.— Schematic of the anisotropic condensation process. The blue region (top left) shows a 
section of the post-shock layer created by a converging flow. Contraction initially begins along the 
direction of the post-shock magnetic held, which is nearly parallel to the post-shock layer. This 
contraction is indicated for a cylinder of initial length L and radius R (top right). When the length 
of the cylinder has shrunk to satisfy V ^ 2i?, it can be treated as an isotropic sphere with radius 
i?, which will collapse if the self-gravity overcomes thermal pressure. Contraction along the post¬ 
shock magnetic held creates dense filaments, and the densest regions within the filament continue 
contracting as quasi-spherical cores (bottom left). 


it is gravitationally stable to transverse contraction across the magnetic held (Mestel & Spitzer 


1956). However, the magnetic held does not prohibit contraction along the length of the cylinder. 


and gravity will be able to overcome pressure forces if L exceeds the thermal Jeans length within 
the post-shock layer, Ly2D = or Lj^sd = ^^(Tr/Gpps)^/^. Here, Dps is the total surface 

density of the post-shock layer. In this situation, longitudinal contraction along a flux tube can 
continue until an approximately isotropic core with L' ^ 2R is produced (Figure]^ bottom)^ with 
density 

, L 


P - 2^Pps. 


( 2 ) 


At this point, transverse contraction is no longer impeded by the magnetic held provided the original 
L ^ Tmag,crit SO the core is magnetically supercritical (note that the mass-to-magnetic flux ratio 
remains the same during the longitudinal contraction). The core will also have sufficient gravity to 
overcome thermal pressure support provided its mass is comparable to that of a critical Bonnor- 
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Ebert sphere at ambient density p', which corresponds to radius (prior to central concentration) 


R ~ -Rth,sph = 2.3 


V AirGp' 


(3) 


Combining L ~ Emag.crit with Equations (Q-Q, this yields 


TD 2 

p' = 0.19 P" 


dvrc. 


0.38 


Pqvq 


(4) 


In Equation Q, we have assumed a strong magnetized isothermal shock with downstream 
magnetic pressure balanced by upstream ram pressure {Bps^/ (Stt) so that 


Bps = 31.04 pG 


Vo 


1 


( 


no \i/2 


km/s ) VlO^cm ^ ) 


(5) 


where po^ vq are the density and inflow velocity of the shock, respectively, and no = Po/l^n foi* 
Pn = 2.3mH the mean molecular weight. We can then solve for the critical radius and mass that 
allows an anisotropically formed core to be both magnetically and thermally supercritical: 


Rc.r\t =5.3 


VGBps 

— 0.09 pc 


= 1.06 


no \-V2 / VQ \ 
/ \1 km/s / 


-1 


1000 cm 


-3 


0.2 km/s 


( 6 ) 


and 


^crit — 14 


G3/2R 


= 2.^ 


ps 


= 2.1 Mrr 


■o 


a/ G^pqVo'^ 
no ^-1/2 Vq 


1000 cm 


-3 


km/s 


-1 


\ 4 

0.2 km/s/ 


(7) 


Equation Q uses Merit = 7ri?crit^Emag,critPps = Rcrit^Bps/{2t/G). Equations ([^ and Q suggest 
that the characteristic mass of prestellar cores formed in post-shock regions in magnetized GMCs 
is determined by the dynamical pressure in the cloud, independent of the cloud’s magnetization, 
when anisotropic condensation along the magnetic field is taken into account. CQ14 already showed 
that models with varying upstream magnetic field directions have similar values of the median core 
mass and radius. Here, we extend our previous investigation to consider variation in the inflow 
velocities and background magnetic field strength. 


^The post-shock magnetic, thermal, and dynamic pressures can be directly measured in the simulation. From 
values listed in Table 1, the magnetic pressure in the post-shock region (Bps^/(87r)) is comparable to the upstream 
ram pressure (povo^), while the thermal pressure (ppsCs^) and dynamic pressure (pps'Crms^) in the post-shock layer 
are about 1 — 2 orders of magnitude smaller. Therefore it is safe to say that the post-shock magnetic pressure is the 
dominant component that balances the upstream ram pressure. 






















3. Numerical Methods and Models 


The simulation setup is similar to C014, and is summarized here. We employ a three- 
dimensional ideal MHD model with convergent flow, self-gravity, and a perturbed turbulent ve¬ 
locity field (Gong & Ostriker 2011). These numerical simulations are conducted using the Athena 
MHD code ( Stone et al.||2008 ) with the Roe Riemann solver. As we found in COM that ambipolar 
diffusion plays a secondary role in core formation, here we consider ideal MHD. The equations we 
solve are: 


| + V.(pv)=0, 

(8a) 

9^^;'+V.(,w-®f)+VP-=0, 

(8b) 

SB 

— + V X (B X v) = 0, 

(8c) 


where P* = P -h P^/(87r). For simplicity, we adopt an isothermal equation of state P = with 
Cs — 0.2 km/s. For both the whole simulation box initially and the inflowing gas subsequently, 
we apply perturbations following a Gaussian random distribution with a Fourier power spectrum 
v‘^ (k) oc (Gong & Ostriker 2011). The amplitude of the velocity dispersion 6v thus depends 
on the Mach number of the inflow M (see Equations (21) and (22) in COM), as 6v = 0.14 km/s • 
(X/10)V2. 


Our simulation box is 1 pc on each side, representing a region within a GMC where a large- 
scale supersonic converging flow with velocity vq — vq z and —vq (i.e. in the center-of-momentum 
frame) collides. The 2 ;-direction is the large-scale inflow direction, and we adopt periodic boundary 
conditions in the x- and p-directions. We initialize the background magnetic field in the cloud, Pq, 
in the x-z plane, with an angle 9 = 20° with respect to the convergent flow. The number density 
of the neutrals, defined as n = p/pn, is set to no = 1000 cm“^ in the initial conditions and in the 
upstream converging flow. The physical parameters defining each model are then the inflow Mach 
number and upstream magnetic field strength M. = v^^/cs and Bq. We choose A4 = 5, 10, and 20 
to look at the dependence of core mass/size on the inflow velocity, and Pq = 5, 10, and 20 pG to 
test whether the initial magnetization of the cloud can affect the core properties (see Table [^. 


Similar to our previous work, we repeat each model 6 times with different random realizations 
of the same turbulent power spectrum to collect sufficient statistical information. Note that the 
resolution adopted in COM| (Ax = 1/256 pc) is not high enough to resolve strong shocks generated 
by high inflow velocity, especially A4 = 20 cases. Therefore, we increased our resolution to 512^ 
for all models in this work, such that Ax ^ 0.002 pc. We use H-correction ( Sanders et al.|[T9M ) to 


suppress the carbuncle instability, and, when needed, first order flux correction (Lemaster & Stone 


2009) (i.e. dropping back to first-order flux estimates for updating the gas variables, if higher-order 


estimates results in negative density). 


2011 ), 


From each simulation, we apply the GRID core-finding method (Gong &; Ostriker 
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Table 1: Summary of the post-shock properties measured and derived from simulations. 



cloud conditions 

simulated post-shock properties* 

corresponding physical scales^ 

Model 

^0 

Bo 

^ps 

Bps 

/dps 

'^^rms 

■^^mag,sph 

-^mag,sph 

Mth.sph 

-^th,sph 

Tmag,crit 


(km/s) 

ipG) 

(10® cm-®) 

ipG) 


(km/s) 

(pc) 

(pc) 

(pc) 

(pc) 

(pc) 

M5B10 

1 

10 

12.1 

43.0 

0.25 

0.44 

7.64 

0.14 

1.26 

0.07 

0.19 

MIOBIO 

2 

10 

18.7 

66.8 

0.16 

0.57 

12.0 

0.14 

1.02 

0.06 

0.18 

M20B10 

4 

10 

32.6 

121 

0.09 

0.79 

23.5 

0.14 

0.77 

0.05 

0.19 

M10B5 

2 

5 

31.8 

59.3 

0.35 

0.41 

2.90 

0.07 

0.78 

0.05 

0.10 

MIOBIO 

2 

10 

18.7 

66.8 

0.16 

0.57 

12.0 

0.14 

1.02 

0.06 

0.18 

M10B20 

2 

20 

9.89 

70.9 

0.08 

0.98 

51.3 

0.28 

1.40 

0.08 

0.37 


*Post-shock properties are measured at t = 0.2 Myr in each model, averaged over the whole post-shock layer. The timescale 
is chosen so the downstream properties are measured before the post-shock layer becomes strongly self-gravitating. 

^See Equations M and (11)-(14) for definitions of physical scales. 


which uses the largest closed gravitational potential contours around single local minimums as core 
boundaries. We then select the gravitationally bound cores as those with negative total energy 
(sum of gravitational, magnetic, and thermal energy). It is then straightforward to measure the 
mass and size for each identified core. For the magnetic flux within a core, we first find the plane 
that includes the minimum of the core’s gravitational potential and is perpendicular to the average 
magnetic field direction within the core. This plane separates the core into an upper half and a 
lower half, and we can measure the magnetic flux $5 through the core by summing up B • fi in 
either the upper or lower half of the core surface (see [COM). The normalized mass-to-magnetic 
flux ratio of the core is therefore F = Mj^B • 27r\/G. 


4. Post-shock Environment and Structure Formation 


4.1. Post-shock Layer 


The post-shock results from our simulations are summarized in Table Similar to COM 


we 


measured the post-shock properties at t = 0.2 Myr, a timescale that is short enough that no cores 
have formed, yet long enough for the post-shock region to reach a steady-state solutionj^ as derived 
in Section 2.1 of C014 In fact, the timescale tgg necessary for the post-shock layer to become 


self-gravitating can be derived by considering when the gravitational weight. 




ttG {2pQVQt, 


sg; 


(9) 


^Although the physical conditions within the post-shock layer do not vary much in time during the initial stages 
of evolution, its thickness grows (approximately linearly in time) because mass is accumulated from the continual 
inflow. The evolution of the layer’s thickness is consistent with the expectation — (2po'^^o/pps)t (see Figure 0. 
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Fig. 2.— The post-shock layer thickness in different models measured from simulations (solid)^ 


compared to possible mass-gathering scales, i?th,sph from Equation (11) (dashed) and i?mag,sph 
from Equation ( [T^ (dotted) within the post-shock layer. Since the post-shock layer is strongly 
magnetized with i?mag,sph rnuch larger than Hpg during the core building phase 0.5 Myr), cores 
cannot collect mass along the direction perpendicular to the layer. 
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exceeds the post-shock pressure Bp^ 


*/87r « PoVq^- The result is 


1 


= 0.79 Myr 


no 


- 1/2 


.1000 cm 

This justifies our choice of measuring post-shock properties at t = 0.2 Myr. 


( 10 ) 


As explained in COM, there are two different length scales (and corresponding characteristic 
masses) for spherical cores in the post-shock region at a given ambient density p: one that is 
supported by thermal pressure (a critical Bonnor-Ebert sphere) 


^th,sph — 2.3 


y/AnGp ’ 


M;h,sph = 4.18- 


y/47rG^p’ 

and one that is supported by magnetic stresses (defined from 4i?/3 = Tmag,crit) 

3 B 


( 11 ) 

( 12 ) 


R 


M 


“ SttVG P ’ 
9 


mag,sph - i287r2G3/2 p 2 


(13) 

(14) 


(see Equations (11)-(12) and (14)-(15) in C014). Figure shows the measured post-shock layer 
thickness in each model, compared with these two possible mass-gathering scales in the post-shock 
environment, i?th,sph and i?mag,sph- H is obvious from Figure]^ that i?mag,sph is much larger than 
the post-shock thickness during the entire core-building phase, and thus magnetically supercritical 
cores cannot form spherically symmetrically within the post-shock layer. Quantitatively, since the 
post-shock layer thickness is = Epg/ (2ppg), we have 


R 


'mag,sph 




ps 


^/87Tp^V^ 


ps 


A7T\/G ^ps 

3 

^/87TGpot 


A7T\/G 

= 1.2 


2povot 


/ no x 

-1/2/ t \ 

VlOOO cm 

yMyiJ 


-1 


(15) 


Since the core formation timescales in our models all satisfy t < 1 Myr, Equation (15) suggests 
i?mag,sph > ^ps when cores formed. This means that gravity-induced mass collection in the direction 
perpendicular to the shocked layer is prevented by magnetic forces, and in-plane mass collection is 
required for core formation in post-shock regions. 


4.2. Structure Formation 


Figure shows examples of structures formed within the post-shock layers, at the time that 
the most evolved core collapses (tcoii; see Section [^. We have selected models with identical 
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Fig. 3.— An example from one of our simulation runs showing the structure formed in the post¬ 
shock layer (in column density; color map) for models with different inffow Mach numbers and 
background magnetic fields. Magnetic field directions in the post-shock layer are also shown {pink 
segments). 
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M5 MIO M20 



X [pc] y [pc] z [pc] X [pc] y [pc] z [pc] x [pc] y [pc] z [pc] 


Fig. 4.— The space-time diagrams of Vx^ Vy^ Vz around the most evolved core in each model, nor- 
malized by the total velocity 'L’tot = ^ each zone. In all models, Vz dominates 

in the beginning of the simulation because of the convergent flow setup, but Vx (along the magnetic 
field lines) soon becomes the strongest component around the forming core. 
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initial turbulence realization, which is responsible for seeding the structures that subsequently 
grow. Filamentary structures are obviously seen in all models with width ^ 0.05 pc, similar to 
those found in observations (see review in Andre et al.||2014 ). Also, note that the filaments are not 
necessary perpendicular to the magnetic field, because the locations of nulls in the velocity field 
are independent of each other on each magnetic field line. 


In addition, we see networks of small sub-filaments or striations parallel to the magnetic field 


in some models. Similar features have been observed in multiple molecular clouds (Goldsmith et 


al. 2008; Sugitani et al. 2011; Hennemann et al. 2012; Andre et al. 2014), and are consistent with 


the theoretical expectation of anisotropy of magnetized turbulence (Goldreich & Sridhar 1995). 


Quantitatively, computational studies suggest /3 < 0.2 is required to have significant anisotropy 


at Mach number = 5 (Vestuto et al. 2003; Heyer et al. 2008), and the critical value of (3 may 


become smaller for higher Mach numbers ( Heyer et al.||2008 ). This roughly agrees with our results 
in Figure striations parallel to the magnetic field direction (not necessarily perpendicular to 
the main filaments) are evident in models with low Mach numbers or strong magnetic fields (M5, 
MIO/BIO, B20). Otherwise, the high velocity turbulence (M20) or the weak magnetization (B5, 
see Table may have destroyed the anisotropy. 

Similar to CQ14[ we use space-time diagrams of different velocity components to demonstrate 
the anisotropic process of core formation (Figure]^. We consider the region with size LxXLyXLz = 
Lmag,crit X 2i?th,sph X 2i?th,sph Centered around the most-evolved core at tcoii of each model, and plot 
the averaged Vx^ Vy^ Vz along x-, ^-, z-directions in the unit of the total velocity 'L’tot- Anisotropic 
gas flows along the x-direction are obvious in all models, and appear much earlier than the core 
collapse (when all three velocity components show convergent flow). Note that, from Figure]^ we 
can see that Model B5 has less prominent convergent flow along the x-direction than the other 
models, indicating that anisotropy is not as strong in this model (see Section]^. 


Quantitatively, if we define overdense (filamentary) structures as those with surface density 
contrast higher than a certain value, say, S > X • Dps, then we can measure the mean surface 
density of filaments, SfUa, as the ratio of total mass inside filamentary structures, 

Mfiia = J S(a;, y) dxdy, (16) 

S> 


to total area (^fiia) of the same structures. The filament formation efficiency (FFE) is defined by: 


FFE = 


Vaia 

Mps 


Mela 
2povot' 


(17) 


Table fists (at time t = tcoll for each model) the measured FFE and total area of filaments 
using X = 1.0 and X = 1.5, as well as three mass-accreting scales Tmag.crit (see Equation 0)> 
Lace, and Xm (see discussion in Section]^. Though the core collapse timescale varies with inflow 
Mach number, the filament formation efficiency and the total area of filaments do not seem to have 
strong dependence on either the inflow Mach number or the pre-shock magnetic field. This is in 
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Table 2: Results from filaments measured at t = tcoii, averaged over all 6 runs for 
each parameter set. 


Model 

^coll^ 

(Myr) 

FFEi.o^ 

FFEi.st 

(pc^) 

7lfila,1.5 

(pc^) 

FmagjCrit 

(pc) 

L ^ 

-^acc 

(pc) 

(pc) 

M5B10 

0.83 

0.65 

0.31 

0.34 

0.11 

0.19 

0.51 

0.39 

MIOBIO 

0.53 

0.57 

0.27 

0.34 

0.11 

0.18 

0.49 

0.28 

M20B10 

0.43 

0.59 

0.31 

0.35 

0.13 

0.19 

0.51 

0.20 

M10B5 

0.58 

0.61 

0.33 

0.34 

0.13 

0.10 

0.27 

0.28 

MIOBIO 

0.53 

0.57 

0.27 

0.34 

0.11 

0.18 

0.49 

0.28 

M10B20 

0.63 

0.53 

0.36 

0.30 

0.15 

0.37 

1.00 

0.28 


^Collapse is defined as the time when nmax = 10^ cm“^ in each simulation. 

^FFE (hlament formation efficiency) is the ratio of the total mass in hlamentary structures 
to the total mass in the shocked layer at tcoii, as dehned in Equation (17). 

^See Section 


contrast to the core formation efficiency (CFE), which varies with tcoii (see Table and discussion 
in Section]^. 

We have also explored how the FFE varies in time (based on individual models). We find that 
the FFE it is fairly constant during self-gravitating stages, with < 15% difference from t — 0.5 tcoii 
to tcoii in all models. This, in addition to supporting our choice of measuring FFE at tcoii, indicates 
that the FFE is not strongly sensitive to the exact age of a cloud. 


Note that there is some arbitrariness in the choice of X. Figure compares the post-shock 
structures under different cutoff values in surface density, and shows the differential PDFs of fila¬ 
ment mass and area as functions of the surface density ratio X = E/Epg. Since there is no “break” 
in the differential PDF at any particular value of X, there is not an obvious value of X to use 
as a lower limit for filament gas. Using X = 1.0, Figure shows X^ia for each model, as well 
as the average column density of the post-shock layer, Xpg = Eps//in, as functions of time. We 
also measured the filament column density at tcoii; in all models, Xcoii = Xfiia(tcoii) ^ 10^^ cm“^, 
comparable to the observed “critical column density” for filaments with active core formation (see 
Andre et al.||2014 ). In detail, we find that for X = 1, N^ow/N^^ 1.8 (see Figure]^. 


review m 


As we shall show below (see Equation (HUi), the expected post-shock column density at the 
collapse time is Xpg oc When E/Epg > X = 1 is used to define filaments, Xcoii/Xpg ps 1.8 

for all models (see Figure [^, implying the same dependence of filament column density on pq ns 
mean post-shock column density, Xcoii oc (see models with different inflow Mach number in 

the left panel of Figure [^. 
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Fig. 5.— Comparison between filamentary structures above different cut-off X values in the crite¬ 
rion S > X • Dps (top), and the fraction of filament mass {bottom left) and the fraction of filament 
area {bottom right) as functions of X, from model MIOBIO. 


5. Statistical Core Properties 


Similar to COM 


we define the timescale at which n^ax ^ lO" 


cm“^ as the moment tcoii 
when the most evolved core collapsesj^ then identify cores formed at this time and investigate 
their physical properties (see Section]^. We note that the cores identified in this way correspond 
to what are termed cores” in |Gong fc Ostriker (2015). Figures and show the statistical 
distributions of core mass, size, mean magnetic field, and mass-to-flux ratio measured from our 
simulations, normalized by total number of cores identified for each parameter set. The normalized 
mass-to-magnetic flux ratio is defined as 


r^X.2nVG. 


(18) 


Cores with F > 1 are magnetically supercritical, and have self-gravity strong enough to overcome 
the magnetic support and collapse. 

Cores identified in our simulations have masses Mcore ^ 0.002 — 10 Mq, sizes i?core 0.004 — 
0.05 pc, and normalized mass-to-flux ratio F ^ 0.4 —4.5, consistent with observations (e.g. Troland 


^We have tested using n — 10^ cm~^ and the results are almost the same, since there is little large-scale evolution 
over that time difference. Based on these tests we found that once a core’s maximum density reaches ^ 10^'^ cm~^ 
it is definitely collapsing, and therefore we used Umax > 10^ cm“^ as the criterion of tcoii in the study. 
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time (Myr) 


Fig. 6.— An example from one of our simulation runs showing the average column density of the 
overdense “filament” structures defined (see text) using X = 1 (solid)^ post-shock layer (dashed)^ 
and theoretical value for the post-shock layer (dash-dotted) defined as Xpg = 2noVot. The core 
collapse time is labeled with dotted lines, with corresponding Xcoii = (Aoii)- 
























Table 3: Results from identified cores measured at t 


Model 

# Cores 

Identified'*' 

CFE^ 

(%) 

(Myr) 

'^core 

(10^ cm“^) 

B ^ 

-L ^core 

(pc) 

^coie 

(Mo) 

Bcore 

(/xG) 

F core 

Mcore/M be 

M5B10 

34 

6.55 

0.83 

2.7 

0.022 

0.81 

49 

2.3 

2.13 

MIOBIO 

30 

3.65 

0.53 

4.9 

0.014 

0.45 

69 

2.1 

2.46 

M20B10 

28 

0.81 

0.43 

11 

0.009 

0.23 

156 

1.2 

1.17 

M10B5 

46 

1.18 

0.58 

7.7 

0.011 

0.25 

89 

2.1 

0.95 

MIOBIO 

30 

3.65 

0.53 

4.9 

0.014 

0.45 

69 

2.1 

2.46 

M10B20 

59 

3.90 

0.63 

6.9 

0.015 

0.55 

103 

1.7 

2.66 


^Columns (2) —(4) are averages over 6 simulation runs for each parameter set. Columns (5) —(10) are median 
values over all cores for each parameter set (6 simulation runs). 

*We only consider gravitationally bound cores with E^grav + ^thermal + < 0. 

^CFE is the ratio of the total mass in cores to the total mass in the shocked layer at tcoii (see Equation (26) in 

COTI ). 

^Collapse is defined as the time when Umax = 10^ cm“^ in each simulation. The tcoii shown here is the mean 
value over all 6 runs for each parameter set. 

^Rcore is calculated from the total number of zones N within a core, for an equivalent spherical volume: 
Rcore = (3A/'/(47r))^/^ Ax, where Ax = 1/512 pc is the grid size. 


& Crutcher 2008; Sadavoy et al. 2010a; Kirk et al. 2013[ ). We also included the normalized mass 


distribution of starless cores in the Perseus molecular cloud (adopted from Sadavoy et al.||2010a ) in 
Figure]^ as a comparison (see Section]^ for more discussion). The median values of core properties 
are summarized in Table|^ as well as the averaged core formation efficiency (CFE) and core collapse 
time tcoii- In Figure we show that the CFE is positively related to the core collapse time, tcoii- 
This is because more structures in the post-shock region have become nonlinear at later time. In 
fact, we found that for individual models, the CFE can increase by an order of magnitude from 
t — 0.8 t^oii to 


Note that though the mean core density, rlcore? is ^ 10 times larger than the ambient density 
in the post-shock layer, the magnetic field within cores (Score) is not significantly different from 
the post-shock region (see rips and Spg in Table [^. This is additional evidence of anisotropic 
core formation: cores gather material along the magnetic field and become more massive without 
significantly compressing the field and enhancing the magnetic support. 

In the anisotropic condensation model (Section]^, core properties are expected to depend 
on the inflow Mach number. In particular. Equations Q and ^ suggest that Score and Mcore 
should decrease with increasing A^, while varying Bq should not have significant effect on these 
core properties. Furthermore, the core field is expected to be comparable to the post-shock value 
given in Equation ([^, so that it increases with Ai and is insensitive to Bq. Our results in Table 
and Figures!^ and [^generally agree with these theoretical predictions. 
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Fig. 7.— Statistical distribution of core mass {left panel) and size {right panel) for models with 
different inflow Mach numbers {top row) and cloud magnetic fields {bottom row). The bin sizes are 
10 ^-^ Mq and 10^-^ pc, respectively. 


Quantitatively, we plot the median values of core mass, size, and mean magnetic field as well 
as the average core collapse time in Figure [To| as functions of initial Mach number {top row) and 
pre-shock cloud magnetic field {bottom row). We also include theoretical models {dotted lines) with 
Mcore ^ (according to Equation Q), i?core (according to Equation Q), Score M 

(according to Equation ([^), and Soil oc (see Equation (23) below). For each theoretical 

comparison, we adopt the predicted scaling and obtain a best-fit coefficient. All simulated results 
fit the theoretical predictions very well, providing quantitative support for the anisotropic core 
formation model. The fit coefficients we find for radius and mass are Mcore = 4.4 Mq and 

Score = 0.14 pc these are shown in Figure [To| 

The Bonnor-Ebert critical radius and mass for an external pressure Pext are given by Sbe = 


0.485 C 32 and Mbe = 1-2 (G^Pext) 


- 1/2 


If we take Pext ^ and normalize to 
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Fig. 8.— Statistical distribution of core mean magnetic field {left panel) and mass-to-flux ratio 
{right panel) for models with different inflow Mach numbers {top row) and cloud magnetic fields 
{bottom row). The bin sizes are 10^'^ /iG and 0.5, respectively. 


no = 1000 cm Cg = 0.2 km/s as in our simulations, the result is 

\ 3 

M-\ 


-RBE,dyn = 0.196 pC { 
-^BE,dyn = 4.43 Mq { 


1000 cm ^ 
no 


2 km/s 

r(s' '■ 


.1000 cm \y0.2 km/s^ 

Comparing to our fitted core radius and mass expressions, we have 

-Rcore — 0.71 i?BE,dyn) Mcore — 0.99 MBE,dyn- 


(19) 


( 20 ) 


Therefore, our results suggest that bound core properties are well described by critical Bonnor-Ebert 
spheres defined by the dynamical pressure of the environment. This snpports the key conclnsion 
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Fig. 9.— The core formation efficiency (CFE) vs. core collapse timescale. Each point represents 
one model parameter set (M5, MIOBIO, etc.). 



predicted in our anisotropic core formation model]^ 

Equations Q and 0 were derived assuming that the accumulation length along the mag¬ 
netic field is ivmag,crit (Equation Q). If, however, we instead assume an accumulation length 
Lace and follow the same steps as before. Equations 0 and 0 would have an additional factor 

(Tacc/Lmag,crit) 5 l-^- 




^core 


= 0.43 pef---^- 3 

VlOOO cm 


^ 0.2 km/s^ 


( 21 ) 


and 


Mcore = 10.5 M 


o 


( _^^ 

VlOOO cm“^/ 


-1/2 


-1 


, . M~^ -. ( 22 ) 

^ 0.2 km/Sy V-^mag,crit/ 

Comparing to our fits, this implies Lacc/^mag,crit = 2.4 or 3.2 for the mass or radius fit, respectively. 
This suggests that cores actually need to gather material along the magnetic field lines from a length 
scale Lace > Tmag,crit- Since Lmag,crit represents the critical (minimum) length scale for cores to 
be magnetically supercritical, our finding of Lace > Lmag,crit is consistent with the anisotropic core 
formation model. Table includes the value (in pc) of Lace = 2.7 Lniag,crit in each model that 
would be required for the median core mass and radius to match Equations (22) and (21). 


^Note that LBE,dyn and MBE,dyn are respectively factors 0.46 and 0.43 smaller than the radius and mass given in 
Equations ^ and Q. 
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Fig. 10.— Summary of simulated core statistical properties for models with different inflow Mach 
numbers (top row) and cloud magnetic fields {bottom row)^ with theoretical predictions {dotted 
lines). The dashed lines in the core size plots {seeond eolumn) indicate the resolution of our 
simulations; Ax ^ 0.002 pc. 


We also use the best-fit coefficients found in Af-models (Figure 10, top row) to derive the 
predicted values {dotted lines) of core mass, size, magnetic field strength, and collapse time for 
S-models (Figure [To| bottom row). Most of the theoretical predictions are in good agreement with 
the simulation results, except the core mass in the B5 model. This is because the B5 model has very 
strong post-shock density compression but only moderate post-shock magnetic field (see Table [^, 
and supercritical cores may form isotropically. This tendency can also be seen in Figures and 
1^ that the structures formed in the B5 model are more randomly distributed compared to other 
models, the anisotropic gas flow is less prominent, and there is less large-scale structure in the B5 
model. 


Figure 10 shows that the core collapse time follows the relationship tcoii oc Ai very well. 
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as predicted in Equation (29) of Gong &: Ostriker (2011). The best-fit coefficient gives 

tcoll = 1.82 Myr (23) 


If we compare with Equation (29) of Gong & Ostriker (2011) (with no = 1000 cm ^ and = 
0.2 km/s), this would imply a maximum amplification in the post-shock region of ln(5E/5Eo)max = 


2.29. The corresponding length scale of the most-amplified mode (see Equation (30) of Gong h 


Ostriker (2011)) is then 


— 


2 v^ 


1/2 


TT 


2.29 


,- 1/2 


non ( ll'O \~^l^ ( VO \ 
°-^^PHl000 cm-3l (ikm/sj 


- 1/2 


(24) 


In most of our models, Xm > Tmag,crit (see Table [^, which means the most-amplified mode would 
be able to form gravitationally bound cores and collapse. In fact, the amplification ln(5E/5Eo) is 
similar (within 30%) for a range of modes with A up to a factor 4 larger than Xm (see Equation (26) 
of Gong & Ostriker ( |201lD ), so it is not surprising that Lace differs from Xm (see Table [^. 


Using the fitted coefficient of Equation (23) combined with the expectation Loll fhe 

predicted post-shock surface density at the time of collapse is Epg (Loll) = 2po'^o%oin corresponding 
to column density 


^ps (%oii) — 5.4 X 10 


21 -2 
cm 


( _ ^2 _) 

VlOOO cm 


1/2 


^0 


1 km/s 


1/2 


(25) 


This is in good agreement with measured values, as shown in Figure]^ Considering Equation (25) 
and the fact that A/fPa (Loll) /^ps ~ 1-8 in nil models (see Figure|^, this suggests that the filament 
column density at the core collapse time may have the same dependence on inflow density and 
velocity as the post-shock column density, i.e. A^^ia (Loll) (n^o'^o)^^^- 


6. Comparison to the Perseus Molecular Cloud 
6.1. Cloud Environment 

The dark cloud in Perseus is an active star forming region approximately 250 pc away, with a 
total mass of about 10^ Mq over a region about 8 x 25 pc (see review in Bally et al.|2008 ). Dense gas 
tracers and dust emission have revealed filamentary structures and a wealth of dense cores in this 
region (e.g. Enoch et al.||2006 ; Kirk et a'L]|2006D . In addition, since the Perseus molecular cloud has 
been observed in and emission lines (e.g. Ridge et al. 2006), the cloud density should 


be > lO"^ cm similar to the value adopted in our simulations. The Perseus molecular cloud thus 
represents a good case to compare with our simulation results. 
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However, the Perseus molecular cloud shows large velocity differences across the region (Bally 


et al. 2008). The observed CO linewidth is about 5 km/s over the whole cloud (Ridge et al. 
2006). Though numerical simulations with rms Mach number M. — Q — 8 have shown agreement 


with observational data on linewidth and cloud structures ( Padoan et al.|[T9M , 2006), there is still 
uncertainty in the actual value of ay in the Perseus molecular cloud because of the possibility of 


superposition of multiple clouds (Bally et al. 2008). 


For our comparisons, we adopted the observed properties of starless cores in the Perseus 


molecular cloud from Sadavoy et al. (2010a). The core mass distribution of Perseus is included 


in Figure as a comparison to simulations. As discussed in Section the Gaussian-fit peaks of 
the core mass functions from our simulations shift with the inffow Mach number, or equivalently, 
the velocity dispersion in the cloud. From Figure the CMF of Perseus has a peak core mass 
similar to that of the M5 model, suggesting that Perseus may be a relatively quiescent star-forming 
environment with converging ffow velocities only of order ^ 1 km/s. 


6.2. Bonnor-Ebert Mass 


One interesting feature of the Perseus cloud is the existence of “super-Jeans mass cores” (Sa¬ 


davoy et al. 2010b). These massive cores have relatively strong self-gravity compared to their 


internal thermal pressure, but still remain starless. An interesting possibility is that these and 
similar cores may be partially magnetically supported. Our models are useful for addressing this 
question, because we can measure the fraction of super-Jeans mass cores under different environ¬ 
ments in our simulations, and we also can measure magnetic support. 

For consistency with theoretical work, we will consider the critical Bonnor-Ebert mass instead 


of the Jeans mass. We thus convert from the MjMj ratios in Sadavoy et al. (2010b) to M/Mbe, 


making use of the core mass and effective radius published in Sadavoy et al. (2010a), and using 


Equation (19) in Gong & Ostriker (2009): 


c ^ 


Mbe = 1 . 18 — . I - 

C.3 i?3/2 


: 1.18- 

: 1.85 = 

V G Pmean 


1.85 , / : 


G3/2 Ml/2- 


(26) 


For a core at mass M, radius i?, and density Pmean that is pressure confined at its surface, the 
thermal pressure is insufficient to prevent gravitational collapse if M > Mbe- For each core 
identified in our simulations, we calculated the value of the critical Bonnor-Ebert mass using the 
core’s mass and radius. 


Figure 11 shows the statistical distribution of Mcore/MBE from both our simulations and 


Perseus; in addition to the binned counts, we also show best fit lognormal functions for each 
model. The distributions for all models and for Perseus are similar. Figure shows that al- 
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Mcore/ -^BE 


Fig. 11.— Statistical distribution of the ratio between core mass and theoretical Bonnor-Ebert 
mass, compared with the observed values in Perseus, for cores formed in the M-models {top) and 
B-models {bottom). 


though the median core mass is close to Mbe, the majority of our gravitationally-bound cores have 
Mcore/^BE > I 5 naively consistent with the fact that these cores are magnetized. However, these 
super-BE mass cores do not in fact seem to be supported primarily by the magnetic field. Fig¬ 
ure 12 shows the mass-to-flux ratio F versus M/Mbe for all cores from our simulations. Evidently, 
most cores with high M/Mbe 3) are also strongly magnetically supercritical (F > 2 ). This 
suggests that the super-BE mass cores observed in Perseus may be strongly self-gravitating and on 
their way to collapse, rather than being magnetically supported. In fact, in our model, cores with 
Mcore/M be ^ 7 all have nmax ^ 10 ^ cm“^, which means they are the most-evolved collapsing cores 
in individual simulation runs. 
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Fig. 12.— Scatter plot of core mass-to-flux ratio vs. M/Mbe in different models. Each point 
represents one core formed in the corresponding model. 


6.3. Mass-radius Relation 


Several observations have found that there is a power-law relationship between the core mass 
and its size, M oc with k ^ 2.4 (Kirk et al. 2013). Figure 13 is the binned mass-size plot 


from all identified cores in our simulations, compared to the observed cores found in the Perseus 
molecular cloud (reported in Sadavoy et al. 2010a). Similar to the observations, the binned data 
from our simulations show k ^ 2 for the power-law relationship between core mass and radius. At 
a given radius, our cores have slightly higher mass than those in Perseus. 

A relationship Mcore oc i?core^ would suggest that the core surface density Score = ^core/(^ri^core^ 
is constant for cores regardless of their masses and sizes. Figure shows the scatter plot of the 
core column density (Acore = '^core/k^n) versus core mass for all cores formed in our simulations. 
Although core mass varies over nearly three orders of magnitude 0.01 — 10 Mq), Acore is within 
a factor of 10. The mean value is Acore = 3.7 x 10^^ cm“^. By comparison, we found in Section 4^ 
that the overdense filamentary structures have column density AfUa ^ 10^^ 
collapse. Thus, the typical core column density Acore ^ 4 A^ia. 


cm ^ at the time of 


However, any k > 2 value indicates that Acore increases with Mcore or i^core, and this trend is 
evident in Figure [Mj for different models. Figure \lE\ shows the mass-radius relations for individual 
models in our simulations as well as the fitted Mcore Rcore^ power-law (the complete fitting 
coefficients are listed in Table Q. We found that the fitted k values are generally higher than 2, 
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0.001 0.010 0.100 
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Fig. 13.— The mass-radius relationship measured from our simulations (diamonds) compared with 
the observation results from Perseus (asterisks)^ using the median values of the binned counts. For 
both the simulations and observations, the vertical bars represent the ±25% values in each bin. The 
best-fit power laws (dotted lines) are M oc for Perseus, and M oc R?'^^ for our simulations. 


implying that A^core is not a constant over cores with different masses and sizesj^ This means that 
it is possible that there is no “universal” core column density, but simply a weak dependence of 
A^core on parameters, which is difficult to identify from the present models. For example, the post¬ 
shock column density at the time of core collapse varies as A^pg oc (see Equation ([^), 

and filament column densities appear to follow a similar trend. If the mean core column density 
is also a multiple of this, then it would vary by only a factor two for our models, which all have 
no = 1000 cm“^ and have V{) varying by a factor four. We do indeed find a higher mean A^core for 
Vo — A km/s (4.8 x 10^^ cm“^) compared to 'L’q = 1 km/s (3.0 x 10^^ cm“^). Further investigations, 
both observational and computational, are needed to reach a clearer conclusion. 


^Composite distribution of cores from different models show a smaller value of /c, and more dispersion, than 
individual models. 
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Fig. 14.— Scatter plot of A^core vs. Mcore in different models. Each point represents one core 
formed in the corresponding model. 


7. Summary 


In this paper, we extended the investigation of C014 to further examine the anisotropic core 
formation model and test the theoretical scalings of core properties over a larger parameter space. 
We carried out fully three-dimensional ideal MHD simulations with self-gravitating gas, including 
supersonic convergent flows with local turbulence. Our models allow for varying inflow Mach 
number and magnetic field strength of the background cloud. Our simulation results demonstrate 
that the ram pressure of the converging flow (po'^o^) is the dominant factor controlling the physical 
properties of cores formed in the shocked layer. These core properties are consistent with the 
predictions of the anisotropic core formation theory. Although the post-shock layer is strongly 
magnetized in all cases, core properties are insensitive to the pre-shock magnetic field strength. 
We also compared cores formed in our simulations with those observed in the Perseus molecular 
cloud, and found very similar core mass distribution, super-Bonnor-Ebert mass ratio, and mass-size 
relation. 


Our main conclusions are as follows: 


1. Considering typical CMC conditions, spherically symmetric core formation is impossible in 
the magnetized post-shock region, because the required mass gathering scales are much larger 
than the thickness of the shocked layer (Table and Figure [^. Quantitatively, it takes 
> 1 Myr for the post-shock layer thickness to be comparable with the magnetic critical 
length under post-shock conditions (Equation ([T^), much longer than typical core formation 
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Table 4: The fitted Mass-size relationship from 
our simulations, Mcore/^O = ^(^core/pc)^. 


Model 

A (xlO^) 

k 

R-square 

M5 

33.57 

2.83 ±0.18 

0.94 

MIO 

10.28 

2.40 ±0.23 

0.82 

M20 

36.14 

2.58 ±0.21 

0.86 

B5 

4.31 

2.17 ±0.22 

0.64 

B20 

17.66 

2.47 ±0.12 

0.88 

All 

4.45 

2.17 ±0.08 

0.81 


timescale in our simulations. 


2. Filamentary structures formed in the post-shock regions are similar to those found in obser¬ 
vations, with dense cores embedded within filaments (Figure]^. We measured the filament 
formation efficiency (FFE) to be around 50% (dependent on the choice of column density 
threshold of filament; Table [^, independent of the pre-shock conditions. We also found that 
the filament column density at the time when cores start to collapse is proportional to the 
mean post-shock column density; A^ma (Wi) ~ l-8A^ps (Figure]^. 

3. Our velocity space-time diagrams (Figure show clear evidence that the mass-gathering 
flows that create cores and filaments are highly anisotropic. Until late times, flow along 
the magnetic field is much stronger than in the two perpendicular directions. However, our 
simulations also show that the “seeds” of cores are present even at early times. This suggests 
that core and filament formation is simultaneous, instead of the commonly-assumed picture 
that cores form only after filaments do. 


4. Magnetically supercritical cores form within the post-shock layers in all of our simulations, 
with masses ^ 0.002 — 10 Mq, sizes ^ 0.004 — 0.05 pc, and normalized mass-to-flux ratio 
^ 0.4 — 4.5 (Table 1^. The core formation timescale is tcoii ^ 0.4 — 0.9 Myr, and the core 
formation efficiency is positively-related to the core collapse time (Figure |^. 


5. The statistical distributions of core mass, size, mean magnetic field, and mass-to-flux ratio 
clearly show that median core properties depend on the pre-shock inflow Mach number A4 = 
'L’o/cs but not the upstream magnetic field strength Hq (Figures and |^. The theoretical 
scalings predicted in the anisotropic core formation model are Mcore A4“^, i?core A4“^, 
and Score ^ Sps oc M (Equations @-0). which agree with our simulation results very well 
(Figure [To|). Furthermore, the core collapse timescale in our MHD simulations generally follow 


the relationship Soil oc A4 The Soil scaling is consistent with the prediction of Gong fc 


Ostriker (2011) based on hydrodynamic analysis, because the flows in the post-shock layer 
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are primarily parallel to the magnetic field. This also gives the post-shock column density at 


tcoii to be A^ps (tcoii) (Equation (25)). 


6. Quantitatively, the median core mass and radius depend on inflow velocity as Mcore = 
0.88 Mq('?;o/ (km/s))~^ and i?core = 0.028 pc('?;o/ (km/s))~^. This suggests that the core 
mass and radius will be, respectively, a factor 0.99 and 0.71 lower than the Bonnor-Ebert 
critical mass and radius computed using the sound speed and total dynamical pressure (po'^o^) 


in the cloud (Equations (19) and (20)). This result is similar to the scaling for characteristic 


mass proposed by [Padoan et al. (1997), but our measured coefficient is higher by a factor 


~ 2 . 

7. Cores identified in our simulations have physical properties very similar to those observed 
in Perseus ( [Sadavoy et al. 2010a). In addition, we found similar statistical distributions of 
Mcoyc/Mbe in simulations and observations (Figure [TT|). We suggest that the “super-Bonnor- 
Ebert mass cores” identified in Sadavoy et al.| ( |2010bD are probably not supported by magnetic 
pressure and will collapse gravitationally, since most cores with high Mcore /^be in our sim¬ 
ulations also have high T values, indicating that these cores are magnetically supercritical 
(Figure . 

8. We find (Figure [T^ a composite mass-radius relation for our simulated prestellar cores com¬ 
parable to that seen in observations, Mcore Rcore^ with k = 2 — 2.5 (e.g. Kirk et al. ]|2013| ). 
Although the observed relation is sometimes interpreted as implying a “universal” core sur¬ 
face density, our results suggest that there might be a weak dependence of the core surface 
density Score = Mcore/ (j^Rcore^) on corc mass or radius (Figure [l4|). We also find that the 
exponent k in the mass-size relation M oc is larger for individual models with consistent 
shock conditions than the composite from heterogeneous environments (Figure [T^. 


To conclude, the success of the anisotropic core formation model for explaining idealized con¬ 
verging turbulent magnetized flows is very encouraging, and provides strong motivation for testing 
these ideas in global MHD simulations of star-forming molecular clouds. Further investigations 
considering more extreme conditions of GMCs would also be interesting to examine the properties 
of core-forming filaments, and potential variations in the core mass-size relationship. 
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